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ABSTRACT 


We have calculated the O(&ew M? /M2,) supersymmetric electroweak corrections to 
single top quark production via qq’ — tb at the Fermilab Tevatron in the minimal 
supersymmetric model. The supersymmetric electroweak corrections to the cross 
section are a few percent for tan@ > 1, and can exceed 10% for tan@ < 1. The 
combined effects of SUSY electroweak corrections and the Yukawa corrections can 
exceed 10% for favorable parameter values, which might be observable at a high- 


luminosity Tevatron. 
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1. Introduction 


The top quark has now been discovered by the CDF and DO collaborations at the Fer- 
milab Tevatron[1]. Measurements of its mass are 176 + 9GeV and 170 + 18 from CDF and 
DO, respectively, and the world-average value of the top quark mass from Run I at the 
Tevatron was recently reported to be 175 + 6GeV, based on roughly 100 pb~! of data[2]. 
At the Tevatron the dominant production mechanism of the top quark is the QCD pair 
production process gg — tt|3]. However, because the top quark is so heavy, electroweak 
production of single top quarks can become significant, particularly at the next Tevatron 
Run-II. With ys = 2 TeV and an integrated luminosity of 2fb~' one can expect, in the 
Standard Model(SM), that for a 175GeV top quark, there will be about 1.4 x 10ftt pairs 
and 5 x 10° single top events produced|4], which is about 35% of the total tt rate. After 
taking into account the b-tagging efficiency and the detection efficiency[5], there are about 
1000 single-b-tagged tt pairs in the /+jets sample, 100 in the dilepton sample, and 250 single 
top events in the /+jets sample available for testing various properties of the top quark. 
Even with fewer events, single top production processes are important because they involve 
the electroweak interaction and, therefore, can probe the electroweak sector of the theory, 
in contrast to the QCD pair production mechanism, and provide a consistency check on the 
measured parameters of the top quark in the QCD pair production. At the Tevatron sin- 
gle top quarks are produced primarily via the W-gluon fusion process|6] and the Drell-Yan 
type single top process, gg’ — W* — tb(W* process)[7], which can reliably be predicted in 
the SM. The theoretical uncertainty in the cross section is only about a few percent due 
to QCD corections|8]. As analysed in Ref.[9], the statistical error in the measured cross 
section for the W* process at the Tevatron will be about +30%; however, a high-luminosity 
Tevatron would allow a measurement of the cross section with a statistical uncertainty of 
about 6%[9]. At this level of experimental accuracy a calculation of the radiative corrections 
in the SM is necessary and effects beyond the SM, for example SUSY corrections, should 
also be considered. 

In Ref.[9] the QCD and Yukawa corrections to the W* process have been calculated 
in the SM. In a previous paper [10] we calculated the Yukawa corrections to this process 
from the Higgs sector in the general Two-Higgs-Doublet Model(2HDM) and the Minimal 
Supersymmetric Model(MSSM)|11] and found that the corrections can amount to more than 
a 15% reduction in the production cross section relative to the tree leve result in the 2HDM, 
and a 10% enhancement in the MSSM. However, in the MSSM, in addition to these Yukawa 


corrections from the Higgs sector, the Supersymmetric (SUSY) corrections due to super 
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particles (sparticles) should also be taken into account. The dominant virtual effects of 
sparticles arise from the SUSY electroweak corrections of order O(Q¢,M?/Mj,) and the 
SUSY QCD corrections of order O(a) which arise from loops of charginos, neutralinos and 
squarks, and gluinos and squarks. It is well-known|[12] that the anomalous magnetic moment 
for a spin 1/2 fermion vanishes in the SUSY limit and away from the SUSY limit there is a 
partial cancellation. Therefore, in general, one can expect the Yukawa corrections from the 
Higgs sector and the SUSY electroweak corrections from virtual charginos and neutralinos 
to cancel to some extent. 

In this paper we present the calculation of the O(a,M7?/Mj,) SUSY electroweak cor- 
rections at the Fermilab Tevatron in the MSSM, and show the combined effect of including 
all the O(@..?/M;,) terms from both Yukawa corrections and SUSY electroweak correc- 
tions. The O(a,) SUSY QCD corrections will be given elsewhere[13]. The paper is organized 
as follows: In Sec. II we present the analytic results in terms of the well-known standard 
notation of one-loop Feynman integrals. In Sec. III] we give some numerical examples and 


discuss the implications of our reults. 


2. Calculations 


2.1 Formalism 


In our calculations we used dimensional regularization to control all the ultraviolet di- 
vergences in the virtual loop corrections and we adopted the on-mass-shell renormalization 
scheme[14]. Including the O(a¢,,M7?/Mj,) electroweak corrections, the renormalized ampli- 


tude for gq’ — tb can be written as 
Mren = Mo J 5 MSUSY-EW (1) 


where Mo is the tree-level matrix element and 8M USY-EW represents the SUSY electroweak 
corrections. The tree-level Feynman diagram for single top quark production via qq’ — tb is 
shown in Fig.1(a). The amplitude Mo is given by 
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.g E as E 
Mo = i> 0(P2) y Piu(pi)U(p3)y"Prv(pa)- (2) 


The amplitude /M°USY-EW is given in the following section. Here p; and pz denote the 
momentum of the incoming quarks q and q’, while p and p4 are used for the outgoing t and 
b quarks, and ŝ is the center-of-mass energy of the subprocess. 
The renormalized differential cross section for the subprocess is 
da 8 — M? 
T aag Dal Mrenl”, (3) 
where @ is the angle between the top quark and incoming quark. Integrating this differential 


cross section over cos@ one obtains the cross section for subprocess 
ô = ĉo + Ac (4) 
where the tree-level cross section ĉo is given by 


fy = <2 8 Meg gy? + (6 — M2)(M2 MMMA (5) 
T287 26 — M22 3 


A 


and AG represents the SUSY electroweak corrections. 
The total hadronic cross section for the single production of top quarks via qq’ can be 


written in the form 


o(s) = 5 | deidr26,j(x028, MÈ, [f (wr, u) fP (@2, 0) + (A > B)], (6) 
i,j 

where 
s = (A+ By, (7) 
§ = 2X 1X28, (8) 
p = tħ, (9) 

and 

P2 = LP». (10) 


Here A and B denote the incident hadrons and P, and P, are their four-momenta, while 
i,j are the initial partons and x and gz, are their longitudinal momentum fractions. The 


functions ff and fr are the usual parton distributions|15,16]. Finally, introducing the 


convenient variable T = zı% and changing independent variables, the total cross section 


becomes 
ol) =D f CEEE) (11) 


ap 7 T 8 aT 


where 7) = (M; + M,)?/s. The quantity dL;;/dr is the parton luminosity, which is defined 
to be 


Se = [SUA mPa) + (A B) (12) 


2.2 SUSY electroweak corrections 


The SUSY electroweak corrections of order &ewM?/Mê2, to the process qq’ — tb arise 
from the Feynman diagrams shown in Fig.1(b)-(g). The matrix element for these corrections 
can be written as 


5 MSUSY-EW = 5 MSUsy— EW + 6Mbow + 6M, 


- (13) 


bos 


where 6 MSUSY-EW represents corrections arising from the self-energy diagrams and vertex 


Web 
diagrams [Fig.1(b)-(e)], while 6 Myo. and ô MẸ, correspond to the box diagram |Fig1.(f)] and 
crossed box diagram|Fig1.(g)|, respectively. Oy EW is given by 


2 1 1 1 
MPW =i Larl Pa) u Pil Pa Ps) PGZ + 5628 + EL) 
8 — My 


+P} P, By + PP, Es] (P) (14) 


The renormalization constants and form factors in Eq.(14) are 


SZ = aa Pis Pog tE) + 2M Mega Rago 

+MPRaal? + [Laser + 3 sra vaata? L (15) 
e E (16) 

Er = va aA VaR RisgOR Myy Mygco + La gOK Myy Milco + ex) 
Te 9 Kj [ME (c22 — c23) + (c12 + c23) 4 Zea = 31} (17) 
Pe S E a — EAV AROR [LizoMgocia + Rigo Mi(ci2 + c22)), (18) 

and 

Es = a AVAR; |t Or; M xt (co + ex) — Ri gOKj Mile + ¢2s)| (19) 


where sums over i, j, k are implied and the functions c;;(P1, P3, M+, Mg, Myo ) are the three- 
T 


point Feynman integrals|17]. The functions ae i) Got GIK) and constants in the above equations 


are defined as 


y 1 1- k 

POD i dyy? log = uu =A ua y) + mky (20) 
u 
j dui 
QCIk) = -f[ d y ( ; 21 
i o P mbyly — 1) + mad — wy) + my os 
ie. ates Se 22 
Mw sin 8’ (22) 
Ly, = cosĝa Lg = —sind,;, (23) 
Ra = Sines Rọ cos 0s, (24) 
Lag = A; cos 03 — C; sin 64, (25) 
Lins = —A;sin 0; — C} cos 0z, (26) 
Rag = —Aj sin 63 + B} cos 0z, (27) 
Rex = —Aj cos 6; — B;sin 0z, (28) 
1 
Lb + * 
O; = yg alia + NioeV51, (29) 
and i 
Oğ = a Nialie + NigU ya, (30) 
where 6, is the mixing angle of squark q, and 
NF gN! 
A; = E B; =C} + J2 (31) 
2mw sin 8’ 2Cw 
A S EN 

Cy = 3 eN; — 3 Cw Ni (32) 
Ni = NyaCw + NjoSw, (33) 


In the above Sw = sin 0w and Cw = cos @y. The chargino masses M; and matrix elements 
V,; depend on parameters M, u and tan 8, whose expressions can be found in Ref.[11]. The 
neutralino masses Myo and matrix elements N,; are obtained by diagonalising the matrix 
Y [11]. Given the values of the parameters M, M’, u and tan 8, the matrix N and Myo 
can be obtained numerically. Here, u is the coefficient of the Hı — Hə mixing term in the 


superpotential and M and M’ are the masses of gauginos corresponding to SU(2) and U(1), 


respectively. With the grand unification assumption, i.e. SU(2) x U(1) is embedded in some 
grand unified theory, we have the additional relation M’ = 545M f 


The box diagram amplitude 6M,,, is 


Mos = (Ps) |(F Pr + f2Pi+ f? PaPr + ft PuPr) ulPi) 0(P2) PsPr 
+ (FEA Pr + FEY Pr) u(P1) (P2) Pr] v(Pa) (35) 


Here the form factors f?3 5 are 
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f = -JaN VaV pRa LaL [Lig M: (Di2 — Dis + Dea — Drs) 
+ RieMg(Da~ Di) ao) 
2 
f = A Va Vj Re, La Ligo Luze (Dia — Dis + Dos — Das), (37) 
2 
g x x 
f = -gar VAV jR Da Lazgi Par, (a 


and B 4g Can be obtained through the permutation 


b o = f? 39 
F246 fiss Lrg Ry, 30 hg x0 Ra ( ) 


The sums over i, j, k, l are implied and the functions Dj;(P4, P3, — P1, M+, Mi,, Mg, Mg) are 
4 
the four-point Feynman integrals|17]. 


The amplitude for the crossed box diagram Mg, is 
Mios = —tU(Ps) LfrPr+ foPrt+ fs PaPr + fi PoP] u(P2) o(Pi)Pro(Pa), — (40) 


where the form factors f$ are 
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LE -Jra UI Re La Lg Met [MRi (Do + Diz) + MypLizxgDo], (40) 

E = NUV aR La Rag My [MiLa Dot Da) + MaRP], 2) 

fs = -NUV Ri ba Lyg May Rig (Do + Du); (43) 
and 

p= = AU} Vio Re, La Ry ge M+ Ligg l(Do + Du). (44) 


The sums over i, j,k, l are again implied and the functions D;;(P,, P3, —P2, My+, Mi,, Myo, My) 
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are the four-point Feynman integrals[17]. The constants x? Ryze are defined by 
Lax = A; cos0g — C; sin 0z, (45) 
a% = —Ajsin 03 — C; cos 0g, (46) 
Rax = -A7 sin 05 + B; cos 6,, (47) 
and 
ET = — A? cos 05 B; sin 6;, (48) 
with 
Was ging N; puota INja (49) 
2mw cosĝ8 7 1 20w 
and 2 
0; = -eNi + aa N's (50) 


3. Numerical results and conclusion 


In the following we present numerical results for the corrections to the total cross section 
for single top quark production via qq — tb at the Fermilab Tevatron with \/s = 2 TeV. In 
our numerical calculations we used the MRSG parton distribution functions|16] and chose the 
scale u = V8. Also we neglected SUSY corrections to the parton distribution functions. For 
the parameters involved, we chose Mz = 91.188GeV, My = 80.33GeV, M; = 175GeV, Mp = 
5GeV and aew = 1/128. Other parameters were determined as follows: 

(i) The upper bound on tan 8; viz, tan 8 < 0.52GeV~'My+, was determined from data 
on B > rvX|18]. The lower limits on tan 8 are tan 8 > 0.6 from perturbative bounds [19] 
and tan 8 > 0.25 (for Mı = 175GeV) from perturbative unitarity|19]. We limited the value 
of tan Ø to be in the range of 0.25 to 5, as larger values of tan 8 are not interesting, although 


allowed by the current data[18], since the effects are negligibly small. 


(ii) For the parameters M;,, Mz, tan 8 and Mrr = A+ cot 8 in top squark (stop) mass 
matrix[20] 
ij M? +m; + 0.35 cos(28) MZ —m(A; + pecot 8) 7 
ae —m,(A; + u cot 8) MP +m? + 0.16 cos(23)Mz J’ (51) 


we assumed M;,, = Mz,- There are then three free parameters in the stop sector and we chose 
the mass of the lighter stop m; , Mrr and tan 8 to be the three independent parameters. 
The best current lower bound on the stop mass is 55GeV coming from LEP, operating at 
V/s = 130 — 140 GeV[21]. We conservatively took the lower bound to be 50 GeV for m;,. For 
the other squarks; i.e., 7,’ and b, we neglected the mixing between left- and right-handed 
states and assumed Ma = Mg My = Mz = M; Mz, which was then determined 
L 2 1 2 
by[20] 


mg. = m? + Mz. + cos(28)(—5 + 55h) M3 (52) 

(iii) For the parameters M, M', u and tan ĝ in the chargino and neutralino matrix, we 
put M = 200 GeV, u = —100 GeV and then used the relation M’ = 545M [11] to determine 
M'. 

Some typical numerical calculations of the SUSY electroweak corrections are given in 
Figs.2-4. 

Fig.2 shows the SUSY electroweak correction Ac /co as a function of lighter stop mass 
M;,, assuming tan? = 1 and Mrr = m. The correction is only a few percent of the tree- 
level value o and is quite sensitive to the lighter stop mass. There are two peaks at about 
M;, = 75 GeV and 67GeV due to the fact that m; = 175 GeV, Myo = (100, 107, 128, 221) 
GeV and the threshold for open top decay into a neutralino and a lighter stop is crossed in 
these regions. 

Fig.3 gives the SUSY electroweak correction as a function of M;pr, assuming tan 8 = 1 
and Mz, = 60 GeV. This correction is also sensitive to Mrr. With increasing Myr the mass 
splitting between the two stop quarks increases. Since we fixed the mass of ¢,, the mass 
of t then increases with Mrr. Furthermore, with an increase of Myr, Mz, increases and 
thus the sbottom masses also increase, as seen from Eq.(52). Since we assumed the masses 
of the squarks ĝi, q@2,q, and qg5 are degenerate with the sbottoms, these masses then also 
increase with Myr. So, with an increase of Myr, all the squark masses except the lighter 
stop increase and their virtual effects decrease due to decoupling effects. From Fig.3 one sees 
that for Mrr > 200 GeV the magnitude of the correction drops below one percent. 

In Fig.4 we present both the SUSY electroweak correction and the Yukawa correction|10] 
as a function of tan 8. For the SUSY electroweak correction we assumed Mrr = mų. Since 
both these corrections are proportional to we they can be very large for small tan 8. 
From Fig.4 one sees that the SUSY electroweak correction exceeds -10% for tan < 0.5, 
and tends to infinity when tan 8 tends to zero. As in the case of the Yukawa corrections, 


the SUSY electroweak corrections are negligibly small for tan 8 > 1. Also, comparing the 
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SUSY electroweak correction with the Yukawa correction in Fig.4, one notes that the SUSY 
electroweak correction and Yukawa correction have opposite signs, and thus cancel to some 
extent. If the lighter stop quark has the same mass as the charged Higgs boson, the cancel- 
lation is appreciable. However, as seen in Fig.4, if the charged Higgs is much heavier than 
the lighter stop quark, the magnitude of Yukawa correction is much smaller than the SUSY 
electroweak correction and there is very little cancellation. In such a case the combined 
effects can exceed -10% for tan 8 < 1. 

To summarize, the combined effects of SUSY electroweak corrections and the Yukawa 
corrections can exceed 10% for favorable values of the parameters. Since the cross section 
for single top production can be reliably predicted in the SM [9] and the statistical error 
in the measurement of the cross section will be about 6% at a high-luminosity Tevatron|9], 
these corrections may be observable; at the least, interesting new constraints on these models 
can be established. 
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Figure Captions 


Fig.l Feynman diagrams for the SUSY electroweak (SUSY EW) corrections: (a) tree- 
level, (b)-(h) one-loop corrections. 

Fig.2 The SUSY electroweak (SUSY EW) correction Ag/ao as a function of M;,, assum- 
ing tan 8 = 1 and MrR = m4. 

Fig.3 The SUSY electroweak (SUSY EW) correction as a function of Mrr, assuming 
tan 8 = 1 and M; = 60 GeV. 

Fig.4 The SUSY electroweak (SUSY EW) correction and the Yukawa correction as a 


function of tan 8. Mrr = M, was assumed for the SUSY EW correction. 
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